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ABSTRACT 

We present a deep X-ray observation of the young Galactic supernova remnant Gas A, 
acquired with the ROSAT High Resolution Imager. This high dynamic range (232 ks) 
image reveals low-surface-brightness X-ray structure, which appears qualitatively similar 
to corresponding radio features. We consider the correlation between the X-ray and 
radio morphologies and its physical implications. After correcting for the inhomogeneous 
absorption across the remnant, we performed a point by point (4" resolution) surface 
brightness comparison between the X-ray and radio images. We find a strong (r = 0.75) 
log-log correlation, implying an overall relationship of logCE^^^^^) oc (2.21 ± 0.05) x 
log(S^^^.^). This is consistent with proportionate partition (and possibly equipartition) 
between the local magnetic field and the hot gas — implying that Gas A's plasma is fully 
turbulent and continuously amplifying the magnetic field. 

Subject headings: supernova remnants - supernovae: individual (Gas A) — X-rays: ISM 
— turbulence — MHD — techniques: image processing 



The University of Minnesota; E-mail: jonathan@lheamail.gsfc.nasa.gov 
Universities Space Research Association; E-mail: gotthelf@gsfc.nasa.gov 
E-mail: petre@lheavx.gsfc.nasa.gov 



1. Introduction 

A comparison of the X-ray and radio emission of 
young supernova remnants (SNRs) provides a pow- 
erful tool for investigating the physical relationships 
among the thermal plasma, cosmic ray electrons and 
the magnetic field (e.g., Dickel et al. 1982; Matsui et 
al. 1984; Arendt et al. 1990; Dyer & Reynolds 1998). 
Radio emission is governed primarily by the density of 
relativistic electrons and the magnetic field strength, 
while the intensity in soft X-rays is dictated by the 
gas density. Shell SNRs provide excellent laboratories 
for investigating the interaction among these physi- 
cal processes, with Cas A being the natural launch- 
ing point for such investigations, because of its high 
surface brightness in both the radio and X-ray wave- 
length regimes. 

Previous work has shown that the soft X-ray mor- 
phology of Cas A, on angular scales ^30" (or 0.5 pc 
at a distance of 3.4 kpc), is dominated by absorption 
effects (Keohane, Rudnick, & Anderson 1996, here- 
after KRA). When this absorption is taken into ac- 
count, a higher degree of intrinsic correlation between 
the X-ray and radio images is found. Similarly, Holt 
et al. (1994) found a correlation between the radio 
and the hard X-ray (i?>4.5keV) surface brightness, 
which is not affected by absorption. These correla- 
tions demonstrate that a simple relationship may ex- 
ist among the underlying physical parameters. 

In this Letter we investigate the relationship be- 
tween the radio and X-ray emission on smaller angular 
scales. We use a newly acquired ROSAT High Reso- 
lution Imager (HRI) image of Cas A, with a 40 times 
longer exposure than used in KRA. To account for 
the inhomogeneous column density towards Cas A, we 
"deabs orb" the image using H i and OH absorption 
data (§2^). We perform a point by point (4" res- 



olution) surface brightness comparison between the 
X-ray and radio ima ges and find a statistically sig- 
nificant correlation (§2^). We discuss the physical 
implications of this correlation and its implica- 
tions for future observational and theoretical studies 
of Cas A and other young SNRs (§|). 

2. Analysis 

2.1. The Deep ROSAT Image of Cas A 

Deep ROSAT HRI exposures of the SNR Cas A 
were obtained on December 23, 1995 and June 20, 
1996, with the telescope optical axis pointed toward 



the center of the remnant. The standard processing 
yields a total of 232 ks of acceptable observing time, 
which produces the highest dynamic range (500:1) X- 
ray image of Cas A to date. A final image, shown 
in Fig. 1^, was made by cross-correlating and adding 
images from the two exposures, using the radio image 
as an absolute reference. 

In order to reduce instrumental blurring (e.g. the 
HRI "halo") we applied 5 iterations of the Lucy- 
Richardson restoration technique (Richardson 1972; 
Lucy 1974). As a kernel for the deconvolution, we 
used an on-axis analytical radially symmetric point 
spread function (PSF), which is valid because the 
PSF is approximately constant across the 5' extent 
of Cas A. The resolution of the final image was esti- 
mated to be about 4" (0.07 pc). 

To correct any systematic scaling errors that may 
have been produced by the restoration process, we 
used a 4" FWHM Gaussian smoothed version of the 
raw image to rescale the reconstructed image. This 
required a linear transformation with a scale of 0.92 
and an insignificant zero-offset (0.31 cts/beam). 

This new image of Cas A is the first with enough 
sensitivity to reveal X-ray features, on small angular 
scales {6 ^ 4") and low surface brightness (Sjjni ^ 
0.03ctss^^ C)^^), which have been only previously 
detected in the radio. Many of the most striking 
of these features are outside the bright shell, such 
as the northern and southern knots and the eastern 
jet region. To enhance this outlying structure, we 
show a spatially filtered version of the image (Fig. ||), 
produced by dividing the 4" X-ray image by a 60" 
Gaussian smoothed version of itself. This effectively 
removes large scale structure and enhances regions 
away from the bright ring. 

2.2. The Deabsorbed Image 

The large scale differences between CasA's soft 
X-ray and radio morphology are dominated by fore- 
ground absorption effects (KRA). In order to take 
these effects into account, we used H i and OH absorp- 
tion data (Schwarz, Goss, & Kalberla 1997; Bieging 
& Crutcher 1986), which trace the foreground atomic 
and molecular gas respectively. We translated these 
data into an optical depth map using the empirical re- 
lation (specific to Cas A) for each resolution element: 

(1) 
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Fig. 1.— The 232 ks ROSAT HRI raw image scaled 
linearly. The maximum value is 102-^^)^!^. This 

0.5 pixel 

images is oversampled with a pixel size of 0.5". 




Fig. 2. — The Lucy- Richardson deconvolved image 
with 4" resolution. This grcyscale is scaled linearly 
with a maximum value of 5300 . This image is 

4 beam 

oversampled with a 0.5" pixel size, but has a resolu- 
tion of 4". 




Fig. 3.— The deconvolved ROSAT HRI image di- 
vided by a 60" smoothed version of itself. This figure 
is included in order to show the outlying and small 
scale structure revealed by this deep exposure. 
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Fig. 4. — The deabsorbed ROSAT image scaled lin- 
early. This image is oversampled with a 0.5" pixel 
size, but its resolution is 4". Note that the interven- 
ing column density was mapped in absorption, so this 
image is blanked outside of Gas A's radio-bright disk 
where the H l and OH optical depths are unreliable. 
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where is CasA's optical depth over the ROSAT 
HRI bandpass, and EWb_i and EWq-h are the H i 
and OH equivalent widths respectively. A complete 
description of this method is given in KRA. 

Using this HRI-specific 30" resolution optical depth 
map, we produced the deabsorbed X-ray image of 
Cas A shown in Fig. |^. Absorption structure on scales 
between 30" and 4" could affect our results, but the 
OH data of Bieging & Crutcher suggest that clouds 
of scales 30"-60" dominate the absorption structure. 

2.3. A Comparison with the Radio Image 

The brightest features in the X-ray and radio maps 
are for the most part not correlated. However, when 
viewed on logarithmic scales, the deabsorbed X-ray 
map and the radio map appear quite similar (Fig. |^). 
In order to quantify this similarity, we performed a 
point-by-point comparison between the deabsorbed 
X-ray image and the most recent (epoch 1994) A6 cm 
VLA map (Koralesky & Rudnick 1998). The time 
differences between the radio and X-ray observations 
(~1 year) will have no effect, given that typical proper 
motions in Cas A are of order 0.1" yr^^ (Tuffs 1986). 
Figure ^ shows this correlation for 4" pixels. This re- 
lationship has a correlation coefficient of r=0.75, and 
appears to be real but scatter-dominated. 

We can parameterize the relationship between the 
deabsorbed X-ray surface brightness (Exrayo) ^^'^ 
the radio surface brightness (Sradio) as 



oc S 



radio 



(2) 



To find the value of we performed a least-squares 
analysis between the quantities log(I]xiayD) ^^'^ log(Siadio)- 
We found the bisector between the two standard least- 
squares fits, as is appropriate for a scatter domi- 
nated linear correlation between two physical quan- 
tities (Isobe et al. 1990). At 4" binning we find a 
best- fit 7] — 2.21 with a corresponding variance of 
cr^ = 6.8 X IQ-* = (0.026)2. Monte Carlo simulations 
by Isobe et al. suggest that their formula for find- 
ing the variance in the slope underestimates the true 
uncertainty by about 10%, so we have set our "90% 
confidence" errors (Table 0) to the 1.8cr level. 

We performed two simple tests to show that nei- 
ther the correlation nor the scatter are artifacts of 
our analysis. First, to show that the correlation is 
intrinsic to the data, we rotated one image by 180° 
about CasA's geometrical center before re-correlating 
the data, whereby the correlation coefficient dropped 



from r=0.75 to r=0.41. This implies that the correla- 
tion is due to small scale structure, and is not simply 
a byproduct of CasA's overall circular morphology. 
Second, to show the reality of the scatter, we created 
a simulated X-ray map, by assigning to each pixel 
an X-ray surface brightness that combined the radio 
surface brightness of that pixel scaled by a power law 
of 2.21, plus a random term representing the Poisson 
scatter in the X-ray map. The correlation between the 
radio and the simulated X-ray image showed no sig- 
nificant scatter, implying that the scatter in figure 
due to statistical fluctuations is small compared with 
that observed. Possible phys ical explanations for this 
scatter are discussed in S3.2. 



We performed the correlation analysis described 
above on maps with successively larger bin sizes in or- 
der to investigate the robustness of our analysis and 
investigate how increasing the pixel size affects the 
value of rj (Table |^). As the data are binned, the 
best-fit r] decreases, but its maximum value remains 
at about ry,„ax ~ 2.3. 

To illustrate how smoothing could bias the mea- 
sured values of rj, assume that there exists a value 



7yo>l, such that S 



Xrayn 



radio ' 



Thus the av- 



eraged quantities would be subject to the inequality: 



(Exrayo) = (^^ X ^^.^^^.J >Cx (E.adio)'"" 



(3) 



One possible explanation for the decreasing rj with 
resolution is that this bias affects both the dim and 
bright parts of Cas A to about the same degree; how- 
ever in proportion to their average brightness the dim- 
mer regions are affected to a greater degree, so our 
logarithmic fitting returns a flatter slope rj. Another 
possibility is that the binning of data lowers the dy- 
namic range of the fit, thus skewing r/ towards its un- 
correlated value of unity. A third possible explanation 
is that smoothing increases the importance of path- 
length variations as the governing physical parameter, 
thus also skewing 77 towards unity. Regardless of the 
particular explanation, the trend shown in Table |^ 
suggests that the intrinsic value of rj is between 2.2 
and 2.3. 

2.4. Deviations from the Trend 

Figure is a map of the quantity log{'E,'^^^^^ /T,^ ^^^^) , 
or a measure of deviations from the best-fit slope 77. 
Consistent with an overall physical correlation, this 
map shows very little structure. There are no large 
scale features such as the shell or the western hot 
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Fig. 5. — Logarithmically scaled images of Cas A. We present, from left to right, the deabsorbed X-ray image and 
the A6 cm VLA map. The pixel size is 4", the same as the beam size. The images are scales to correspond to the 
axes in Fig ^ 
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Fig. 6. — Log-log plots correlating the images in Fig|^. 
The left-hand plot is a standard scatter-plot with 3948 
data points. For a majority of points, the random er- 
rors are smaller than the symbols plotted. On the 
right is a plot of the median \og[Y,^^_^^^) in evenly 
spaced bins of log{S^^^.^) (triangles) and visa versa 
(squares). Both data sets are scaled to units of their 
respective overall median values. The same lines are 
plotted on each graph: the dashed lines represent the 
two standard least-squares fits (r7xy=1.7, ?7yx=3.0); 
the dotted lines represent the 90% confidence range 
in the slope; and the solid line is the bisector of the 
two least-squares slopes. 




Fig. 7.— Plot of the quantity \og{Y?-^J^JT.^ 
The greyscale ranges from log(min/mean) = - 
(white) to log(max/mean) = 1.67 (black). Like Fig.^ 
this image is oversampled with a 0.5" pixel size, but 
a resolution of 4". 
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Table 1: The X-ray/radio Regression Parameters 



Resolution 


N 


r 


V 


4" 


3948 


0.75 


2.21 ±0.05 


8" 


876 


0.75 


2.19 ±0.11 


16" 


210 


0.76 


2.14 ±0.17 


32" 


47 


0.75 


1.96 ±0.35 


64" 


10 


0.49 


1.43 ±0.72 



N the number of independent points 
r = the correlation coefficient 

77 = the log-log surface brightness slope (see equa- 
tion 0) with 90% confidence errors 



spot. A few key features are present, however. Some 
are relatively radio bright: the northern knot region; 
a knot near the center; and the outlying structure 
in the south. In addition, there is a relatively X-ray 
bright feature on the eastern shell. 

3. Discussion 

3.1. Proportionate Partition as a Possible Ori- 
gin for the Correlation 

The main objective of comparing the X-ray and 
radio emission in Cas A is to identify the physical rela- 
tionships among the thermal plasma, cosmic ray elec- 
trons and magnetic field. In this section we show that 
our results are consistent with proportionate parti- 
tion between the energy content of the gas and mag- 
netic field. We also develop a physical justification 
for this by postulating a turbulent magnetohydro- 
dynamic (MHD) plasma, which we argue is contin- 
uously amplifying the magnetic field. As both the 
trend and the exceptions are real, this scenario should 
be interpreted as valid "on average," but not neces- 
sarily for each feature in Cas A. 

The ROSAT band spectrum of Cas A is dominated 
by line emission from H and He-like ions of Mg and 
Si and the L-shell ions of Fe (Holt et al. 1994). Thus 
by assuming a constant temperature plasma, we can 
relate the X-ray emissivityQ (e^ „y ) to the gas density 
i'lT'gas) as: 

Sxray ^ '^^aS (4) 

Cas A's radio emission mechanism is believed to be 



*Here the emissivity, £xray' defined in an analogous way to 
the radio emissivity, such that oc I e„ dl, where I is 

■y 1 X rayQ J X ray ' 

the path length. 



synchrotron radiation (Rosenberg 1970), arising from 
relativistic electrons (of density n^ci) interacting with 
the magnetic field {B). We therefore characterize the 
radio emissivity {s^_^^^^ ) as: 



l+a 



(5) 



where a is the radio spectral index [a = 0.77, Baars 
et al. 1977). 

Anderson et al. (1994) concluded that the radio 
brightness of the compact features in Cas A is gov- 
erned primarily by magnetic field amplification. For 
this reason, we believe it reasonable to assume a con- 
stant (or uncorrelated) relativistic electron density. 
Therefore, equation ^ becomes: 



oc B 



l+a 



(6) 



In their MHD simulations, Jun, Norman, & Stone 
(1995) show that Ray leigh- Taylor (R-T) instabihties 
in young SNRs are significantly enhanced by the pres- 
ence of a magnetic field. Kelvin-Helmholtz instabil- 
ities, created by the shear fiow along the R-T fin- 
gers, greatly enhance the magnetic field — especially 
on small scales. In addition, MHD studies of simple 
shear systems show that the magnetic field is ampli- 
fied quickly as the energy cascades to smaller scales 
(Malagoli, Bodo, & Rosner 1996). On the scale where 
the magnetic field becomes dynamically important, 
reconnection and dissipation occur. It is therefore 
plausible that the magnetic field in Cas A is being 
amplified to the point of equipartition with the hot 
gas, locally on the smallest observable scales, so: 



B^ 

8^ ^^"'^ 



(7) 



In aggregate this is approximately true. If one as- 
sumes equipartition between the magnetic field and 
the relativistic plasma, an energy density for each of 
these can be derived simply from the overall radio flux 
density. This overall magnetic energy density (■^) is 
comparable to the thermal energy density (J7gas) es- 
timated from the integrated ASCA X-ray spectrum 
(Allen, G.E. private communication). 

This suggestion of equipartition is apparently in- 
compatible with the simulation of Jun & Norman 
(1996), whose 3-D simulations predicted magnetic 
fields significantly lower than equipartition. How- 
ever, they also found that the higher the resolution 
of the simulation (i.e., the lower effective viscosity). 



6 



the greater the resultant magnetic field. With re- 
gard to the analysis presented here, this distinction is 
moot. If the turbulent cascade is halted on the small- 
est scales by the magnetic pressure, the equipartition 
relation (equation |^) should hold. On the other hand, 
if the turbulent cascade is the primary amplifier of the 
magnetic field, but is halted before the magnetic field 
reaches equipartition, then one would expect the sys- 
tem to still be in proportionate partition, i.e.: 

B'^ CX [/gas OC rigas (8) 

The analysis presented here cannot distinguish be- 
tween equipartition (equation 0) and the looser con- 
straint of proportionate partition (equation^). Simi- 
larly, it does not matter here whether Ugas represents 
the turbulent energy density or the thermal energy 
density, so long as the energy density scales linearly, 
on average, with the gas density (t/gas « rigas). We 
now combine equations |[ ^ and ^ to find a relation- 
ship between the X-ray and radio emissivities: 

This simple relationship agrees remarkably with the 
value of rj found in |2.3| . 

A similar and independent conclusion about tur- 
bulent magnetic field amplification was found by Ko- 
ralesky & Rudnick (1998), who examined the changes 
in Gas A's radio emission as a function of azimuth over 
an 11 year baseline. They found that regions with 
the largest decrease in fractional polarization (corre- 
sponding to a decrease in order of the magnetic field) 
are decreasing the least in flux density. They suggest 
that the same process which is randomizing the mag- 
netic field's orientation is also amplifying it, partially 
compensating for expansion losses. 

This relationship between the radio and X-ray sur- 
face brightness is very different from the predictions of 
most other simple models. For example, in an isobaric 
model, the magnetic pressure and thermal pressure 
add to a constant, so the surface brightnesses would 
be anti-correlated. More realistically, a compression 
scenario where the magnetic field scales inversely with 
area so i? oc Jigas, implies that 77 = 1.7, which is also 
inconsistent with the data. If instead, we make the 
equally simple assumption that the relativistic elec- 
tron density is proportional to the thermal gas density 
(nici OC ngas), then we find i] = 2.0. While we cannot 
completely rule out this possibility, we do not favor it. 
The Larmor radii of the relativistic electrons are very 



small ('^10^ cm for a 1 GeV particle in a 1 mG mag- 
netic field), tightly coupling them to magnetic field 
lines, which in turn arc coupled to the plasma. This 
suggests that rirei oc Tigas may be a good assumption, 
but then one would also expect a relation between the 
magnetic field and density. When a magnetic field de- 
pendence is included, rj is significantly reduced. For 
example, if riioi oc rigas oc B then ry = 0.7, which is 
clearly inconsistent with the correlation. 

3.2. Possible Contributions to the Scatter 

The simple model that fits the surface brightness 
ratio rests on three assumptions: constant relativistic 
electron density, constant X-ray emissivity per den- 
sity squared (Sxray/'^gas); ^^'^ proportionate parti- 
tion. The relative flatness of the ratio map shows that 
deviations from these assumptions on large scales are 
absent. It is plausible, however, that the substantial 
scatter of the ratio for individual cells from the trend 
line represents local violation of one or more of these 
assumptions. The log of the rms deviation about the 
trend hne is ±0.38 (a factor of ~2.5); the log of the 
maximum excursion is ±1.5. 

If the scatter is due entirely to variations of rirei or 
-B^/^gas (proportionate partition), we find the corre- 
sponding rms variation of these is ±0.17 in the log, or 
a factor of ~1.5. The X-ray emissivity is a function 
of both gas temperature and metal abundance. The 
emissivity scales linearly with the abundance of the 
metals producing the dominant lines in the ROSAT 
band, namely Mg and Si; thus local rms variations of 
a factor of ^^2.5 in metallicity can account for the scat- 
ter. The emissivity scales with temperature in a less 
straightforward way, but deviations of log(£xray/'T-gas) 
of ±0.4 can be accounted for by log(kT) variations of 
±0.4. 

Scatter of the degree indicated in any of these pa- 
rameters is physically plausible. For instance, small 
scale gas pressure variations associated with a shock 
encountering a density enhancement can be as large as 
a factor of 6 (Dohm-Palmer & Jones 1996), which in 
turn could produce variations of B^/rigas or £x ray /"gas 
within the observed range. Measurements of any of 
these parameters on small angular scales, at best pose 
a challenge to current and future instruments, and at 
worst may be infeasible to measure (e.g., Uici). 
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4. Conclusion 

We have performed a comparison of CasA's radio 
and X-ray emission to a limiting resolution of 4" (0.07 
pc). The strong correlation between radio and X-ray 
surface brightness can be explained by the scenario 
that Cas A has "on average" a spatially constant rel- 
ativistic electron density and proportionate partition 
on small scales between its thermal and magnetic en- 
ergy densities, as would be expected from a fully tur- 
bulent MHD system. 

These results may have important implications for 
theoretical work. The complex plasma physics of 
SNRs often must be simplified in order for computer 
simulations to run in a cost-effective manner. Un- 
fortunately, most often one assumption is spherical 
or cylindrical symmetry, which is incompatible with 
turbulent flow. It may be possible, instead, to use 
the proportionate partition relation suggested by our 
result, along with some other characteristics of fully 
turbulent systems, to produce a more realistic first- 
order description of young, core-collapse, supernova 
remnants. 

We thank Barron Koralesky for kindly providing 
the A6 cm radio data of Cas A in a digitized form. 
We also wish to thank Glenn Allen and John Dickel 
for insightful conversation. This work was supported 
by NASA. 
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